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Summary _ In this study, we were interested in determining if high titered egg adapted modi- 
fied live infectious bronchitis virus (IBV) vaccines contain spike gene related quasispecies that 
undergo selection in chickens, following vaccination. We sequenced the spike glycoprotein of 
12 IBV vaccines (5 different serotypes from 3 different manufacturers) directly from the vaccine 
vial, then compared that sequence with reisolated viruses from vaccinated and contact-exposed 
birds over time. We found differences in the $1 sequence within the same vaccine serotype 
from different manufacturers, differences in $1 sequence between different vaccine serials 
from the same manufacturer, and intra-vaccine differences or quasispecies. Comparing the 
sequence data of the reisolated viruses with the original vaccine virus, we were able to identify 
in vivo selection of viral subpopulations as well as mutations. To our knowledge, this is the first 
report showing selection of a more fit virus subpopulation as well as mutations associated with 
replication of modified live IBV vaccine viruses in chickens. This information is important for 
our understanding of how attenuated virus vaccines, including potential vaccines against the 
SARS-CoV, can ensure long-term survival of the virus and can lead to changes in pathogenesis 
and emergence of new viral pathogens. This information is also valuable for the development 
of safer modified live coronavirus vaccines. 

Published by Elsevier Ltd. 


Introduction 


Coronaviruses are an important group of enveloped single 
stranded RNA viruses that cause highly contagious respi- 
ratory and enteric diseases in a variety of animals and in 
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humans. Coronaviruses belong to the order Nidovirales, fam- 
ily Coronaviridae, and are separated into 3 groups based on 
genetic and antigenic characteristics [1]. Avian infectious 
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bronchitis virus (IBV) belongs to the group 3 coronaviruses 
along with the other avian coronaviruses, whereas coron- 
aviruses in groups 1 and 2 produce disease in many species 
of mammals and in humans. Confusion around the origin of 
severe acute respiratory syndrome coronavirus (SARS-CoV) 
has lead to a proposed group 4; however, genetic studies 
suggest that SARS-CoV diverged from the group 2 viruses, 
thus a group 2b classification was suggested [2]. 

Avian infectious bronchitis virus has an enormous 
economic impact on commercial poultry causing highly con- 
tagious disease in chickens of all ages. Since IBV was first 
described by Schalk and Hawn in the 1930s, reviewed in 
ref. [3], many serotypes have been identified worldwide. 
There are currently 4 reported structural proteins associ- 
ated with the IBV virion; nucleocapsid, membrane, small 
membrane and spike. The spike glycoprotein on the outside 
of the virus contains epitopes associated with serotype dif- 
ferences and binding of neutralizing antibodies, and it plays 
a role in attachment and entry into the host cell. 

Coronaviruses including IBV have been shown to exist as 
a mixture of genetic mutants within an isolate [4,5]. These 
so called quasispecies are generated through genetic muta- 
tions and recombination events that occur when the viral 
RNA-dependent RNA-polymerase, which lacks proofreading 
capabilities, replicates the viral genome [6]. For SARS-CoV, 
non-synonymous and synonymous substitution rates have 
been estimated to be as high as 3.3 x 10-3 and 4.67 x 10-3 
per site per year, respectively [7]. Evolutionary and mutation 
rates for IBV in the face of vaccination were determined to 
be 2.5 and 1.5% per year, respectively, in the hypervariable 
region of the spike glycoprotein; however, in the absence 
of vaccines, IBV lineages appear to evolve at a much slower 
rate (0.3%) [8]. The group 3 coronaviruses, which exist only 
in avian species were found to have a relatively recent diver- 
gence date (1925) and exponential population growth [2]. In 
general, constant population size is associated with natural 
reservoirs of the virus, which remain largely asymptomatic 
in the host, whereas a change from constant to exponential 
population growth indicates emergence of a relatively new 
virus associated with interspecies transmission and disease 
[2]. 

An actively changing genetic population allows rapid 
adaptation to the host through selection of the most fit viral 
subpopulation. Selection ensures long-term survival of the 
virus at the cell, organism, and host population level, and 
can lead to changes in pathogenesis and emergence of new 
viral pathogens [9]. For economic reasons, vaccine manu- 
facturers attempt to produce high titered (>1 x 10® embryo 
infectious dosesg (ElDsq)/ml) IBV vaccines so they can be 
packaged in 5000 or even 10,000 dose vials. The average 
IBV vaccine dose can vary but is generally approximately 
1 x 10* ElDs9. Since the number of intra-isolate genome sub- 
populations are reported to increase dramatically when the 
virus titer increases [9], we were interested in determin- 
ing if quasispecies of high titered modified live IBV vaccines 
contain specific spike gene mutations that undergo selection 
in chickens, following vaccination. In addition, we wanted 
to determine if molecular changes were solely the result 
of selection or if mutations also occur following vaccina- 
tion. 

To study the dynamics of IBV vaccine genetic adaptation 
to the host, we sequenced the spike glycoprotein of 12 IBV 


vaccines (5 different serotypes from 3 different manufac- 
turers) directly from the vaccine vial. Then, we vaccinated 
chickens with those viruses and examined the spike glyco- 
protein sequence of the vaccines reisolated from vaccinated 
and from contact-exposed birds over time. Comparing the 
sequence data of the reisolated viruses with the original 
vaccine virus, we were able to access in vivo selection of 
subpopulations and genetic mutations. 


Materials and methods 


Vaccines 


Twelve IBV commercial vaccines were obtained from three 
different manufacturers designated A, B and C herein, and 
had a titer of at least 1 x 108 EIlD59/ml. The vaccines used 
in this study from manufacturer A were Ark-DPI (Ark/A), 
GA98 (GA/A) and Mass (Mass/A). Vaccines from manufac- 
turer B were Ark-DPI (Ark/B), Mass (Mass/B) and a mixture 
of Mass and Conn (Conn/B). The vaccines from manufacturer 
C were Ark-DPI (Ark/C), DEO72 (Del/C), GA98 (GA/C), Mass 
(Mass/C), a mixture of Mass and Conn (Conn/C), and MassD 
(MassD/C). The lyopholyzed vaccines were resuspended in 
10ml of diethylpyrocarbonate (DEPC) treated water. The 
DE072 vaccine, stored in liquid nitrogen, was thawed and 
diluted using 1 ml of DEPC treated water. 


Birds and housing 


Specific-pathogen-free (SPF) white leghorn chicks (Merial, 
Gainesville, GA) were housed in positive-pressure Horsfal 
isolation units. Feed and water were provided ad libitum, 
and the birds were examined twice daily. 


Experimental design 


To guard against potential cross-contamination, the 12 vac- 
cines were assigned to 4 different groups ensuring that there 
were no vaccines of the same serotype in each group. Two 
different experiments were conducted with vaccine groups 
1 and 2 in the first experiment and vaccine groups 3 and 4 in 
the second experiment. Vaccines used in experiment 1 were 
for group 1 Del/C, GA/C, Mass/A, and for group 2 Ark/C, 
GA/A, and Mass/C. Vaccines used in experiment 2 were for 
group 3 Ark/B, Conn/B, Mass/B, and for group 4 Ark/A, 
Conn/C, and MassD/C. Each vaccine group was tested in a 
separate filtered air positive pressure room containing the 
isolators. One-week old chicks were randomly divided into 
groups of 10 birds for each vaccine type. Five birds in each 
group were inoculated by eye-drop with 30 yl of vaccine 
(equal to one dose, according to the manufacturers rec- 
ommendation), the other 5 birds served as contact-exposed 
birds. A negative control group of 5 birds was included for 
each group of vaccines tested. 

Tracheal swabs from all birds were collected at 3, 6 and 9 
days post-vaccination (dpv), and at necropsy. Birds in exper- 
iment 1 were necropsied at 14 dpv while birds in experiment 
2 were necropsied at 13 dpv. Swabs from 5 vaccinated birds 
and 5 contact-exposed birds for each of the vaccine types 
were pooled separately in 1 ml PBS (pH 7.4). The swabs were 
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stored at —80°C until they were analyzed for IBV RNA by 
RT-PCR. 


Passage in embryonating eggs 


The last positive tracheal swab by RT-PCR from each vac- 
cinated and contact-exposed group of birds was passaged 
a maximum of ten times in 9—11 day old specific pathogen 
free (SPF) embryonating eggs. Three eggs per group were 
inoculated with 0.1 ml of PBS from the tracheal swabs into 
the chorioallantoic sac (CAS) as described [10]. The eggs 
were incubated at 37°C and allantoic fluid was harvested 
from all eggs with live embryos at 48h post-inoculation, 
pooled and kept at —80°C. Subsequent passage was done by 
inoculating 0.2 ml of allantoic fluid from the previous pas- 
sage into the CAS of 9—11 day old embryonating eggs and 
the allantoic fluid was harvested at 48h post-inoculation as 
above. 


Viral RNA extraction 


Viral RNA was extracted from allantoic fluid, rehydrated or 
diluted vaccine, and tracheal swabs using the High Pure RNA 
Isolation Kit (Roche Diagnostics, Indianapolis, IN) according 
to manufacturer’s protocol. 


RT-PCR amplification and cloning of the S1 gene 


Reverse transcriptase-polymerase chain reaction was car- 
ried out using the Titan One Tube RT-PCR kit (Roche 
Diagnostics, Indianapolis, IN) according to manufacturer’s 
protocol with primers NEWS10LIGO5’ [11] and Degenerate3’ 
[12]. Primers specific for the Conn serotype [13] were used 
to screen for the presence of that virus in swabs collected 
from birds vaccinated with vaccine containing both the Mass 
and Conn types and in swabs from the corresponding contact 
birds. The PCR products were agarose gel purified using the 
QlAquick gel extraction kit (Qiagen, Valencia, CA) according 
to the manufacturer’s protocol. To aid in sequence analysis, 
purified PCR products from Conn/B and Conn/C, which were 
mixtures of Mass and Conn vaccines, and Ark/B were cloned 
into the TOPOXL vector using TOPOXL cloning kit (Invitrogen, 
Carlsbad, CA) according to manufacturer’s protocol. The IBV 
real-time RT-PCR analysis on tracheal swabs was performed 
as previously described [4,14]. 


Sequencing and sequence analysis of the S1 gene 


Viral $1 gene sequence was determined for two different 
production lots for each IBV vaccine examined. The S1 gene 
of reisolated viruses from the first and last RT-PCR positive 
tracheal swabs only was also sequenced for both vacci- 
nated and contact-exposed birds. The $1 gene of reisolated 
virus from Ark/A vaccinated birds was only positive by real 
time RT-PCR. The pooled tracheal swabs collected from 
the 13 dpv vaccinated birds were sequenced. In addition, 
each embryo passaged virus was sequenced to determine 
if genetic changes in the reisolated vaccine viruses would 
revert to the original vaccine sequence. Purified RT-PCR 
products or at least 4 cloned RT-PCR products for each virus 


were sequenced using the BigDye® Terminator v3.1 Cycle 
Sequencing Kit (Applied Biosystems, Foster City, CA) accord- 
ing to manufacturer’s protocol. Sequencing reactions were 
purified using Centri-Sep Columns (Princeton Separations, 
Inc, Adelphia, NJ) or with the Performa DTR Ultra Dye Termi- 
nator removal system (Edge BioSystems, Gaithersburg, MD). 
Sequence reactions were examined at the Molecular Genet- 
ics Instrumentation Facility (University of Georgia, Athens, 
GA). 

The S1 sequences from each virus were compiled using 
the EditSeq and MegAlign program (DNASTAR, Inc. Madi- 
son, WI) and open reading frames were predicted using 
ORF finder at National Center for Biotechnology Information 
(NCBI, http://www.ncbi.nlm.nih.gov/gorf/gorf.html). The 
BlastX program (http: //www.ncbi.nim.nih.gov/BLAST/) was 
used to search GenBank (National Center for Biotechnology 
Information, http://www.ncbi.nlm.nih.gov/) for homolo- 
gous IBV S1 sequences. 


Results 


S1 sequence analysis of vaccine viruses 


All nucleotide sequences reported herein have been submit- 
ted to NCBI, GenBank (http: //www.ncbi.nlm.nih.gov/) and 
the accession numbers are listed in Table 1. To determine 
if mutations or selection of a more fit virus subpopulation 
was responsible for the changes observed when IBV adapts 
to an in vivo environment, we first needed to establish the 
exact sequence of the vaccine viruses directly from the man- 
ufacturer. The S1 gene deduced amino acid sequence for 
different production lots of the Ark/A and Ark/C vaccine 
viruses were the same, whereas Ark/B vaccine produc- 
tion lots had 1 amino acid difference (Pro524Ser) in the 
S1 gene. In addition, sequence differences were observed 
between Arkansas vaccines from the different manufactur- 
ers (Table 2). 

There were 26 amino acid differences in the S1 sequences 
between the two production lots for the Conn/B vaccine 
(Table 3). For the Conn/C vaccine, we observed 5 amino 
acid differences between the different production lots. In 
addition, one of the Conn/C vaccine production lots had a 
deletion of amino acid Asn at position 58. 

There were no differences observed between the S1 
sequences obtained from the two vaccine production lots 
for the Del/C, GA/A, GA/C, Mass/A, Mass/B, and Mass/C. 
There were 13 amino acid sequence differences between the 
two production lots of the MassD/C vaccine virus (Table 4). 
MassD/C vial 1 was closest to Mass/A with only 1 amino acid 
sequence difference between them, whereas MassD/C vial 
2 was closest to Mass/C and Mass/B with 6 and 7 amino acid 
sequence differences, respectively. In addition, there were 
differences in the amino acid sequences between the Mass 
vaccines from the different manufacturers. We observed 
14 and 11 amino acid sequence differences between Mass 
vaccines from manufacturer A and the Mass vaccines from 
manufacturers B and C, respectively. We also found 2 amino 
acid sequence differences between Mass vaccines from man- 
ufacturers B and C. Finally, there were 28 amino acid 
sequence differences and one deletion between the GA vac- 
cines from manufacturers A and C (Table 5). 
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Table 1 NCBI? accession numbers for IBV $1 gene sequences obtained directly from virus in vaccine vials, and reisolated vaccine 


virus from vaccinated and contact-exposed birds 


IBV vaccine virus 


S1 sequence source 


NCBI accession number 


Ark/A 


Ark/B 


Ark/C 


Conn/B 


Conn/C 


Del/C 


GA/A 


GA/C 


Mass/A 


Mass/B 


Mass/C 


MassD/C 


4 NCBI = National Center for Biotechnology Information (http: //www.ncbi.nlm.nih.gov/). 


> dpv = day post-vaccination. 


First vaccine lot 
Second vaccine lot 
13 dpv? vaccinated birds 


First vaccine lot 
Second vaccine lot 
3 dpv vaccinated 

9 dpv vaccinated 

6 dpv contact birds 
13 dpv contact birds 


First vaccine lot 
Second vaccine lot 
6 dpv vaccinated birds 


First vaccine lot 
Second vaccine lot 
6 dpv contact birds 
13 dpv contact birds 


First vaccine lot 
Second vaccine lot 


First vaccine lot 
Second vaccine lot 
Reisolated vaccine virus‘ 


First vaccine lot 
Second vaccine lot 

3 dpv vaccinated birds 
9 dpv vaccinated birds 


First vaccine lot 
Second vaccine lot 
Reisolated vaccine virus 


First vaccine lot 
Second vaccine lot 
Reisolated vaccine virus 


First vaccine lot 
Second vaccine lot 
3 dpv vaccinated 

9 dpv vaccinated 

6 dpv contact birds 
13 dpv contact birds 


First vaccine lot 
Second vaccine lot 
Reisolated vaccine virus 


First vaccine lot 
Second vaccine lot 

3 dpv vaccinated birds 
6 dpv contact birds 


EU283045 
EU283046 
EU283047 


EU283048 
EU283049 
EU283050 
EU283051 
EU283052 
EU283053 


EU283054 
EU283055 
EU283056 


EU283057 
EU283058 
EU283059 
EU283060 


EU283061 
EU283062 


EU283063 
EU283064 
EU283065 


EU283066 
EU283067 
EU283068 
EU283069 


EU283070 
EU283071 
EU283072 


EU283073 
EU283074 
EU283075 


EU283076 
EU283077 
EU283078 
EU283079 
EU283080 
EU283081 


EU283082 
EU283083 
EU283084 


EU283085 
EU283086 
EU283087 
EU283088 


© Reisolated vaccine virus indicates that all viruses reisolated from vaccinated and contact-exposed birds had the same sequence. 
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Table 2. Pairwise comparisons of amino acids between Ark 
vaccines from different manufacturers 
Residue position Ark/A Ark/B Ark/C 
43 Tyr His Tyr 
171 Tyr His Tyr 
aS Ser Ala Ser 
249 Phe Ser Phe 
323 Arg Arg Tyr 
410 Lys Glu Lys 
524 Ser Pro (Ser)? Ser 


4 Amino acid in brackets indicates the residue observed in 
vaccine vial 2. 


Virus detection from vaccinated and 
contact-exposed birds 


Tracheal swabs collected from vaccinated and contact- 
exposed birds were tested for IBV by RT-PCR and real-time 
RT-PCR and the results are presented in Table 6. Tracheal 


swabs collected from all the negative control birds were 
negative for IBV by RT-PCR and real-time RT-PCR. Based on 
sequence data from reisolated viruses, there was no cross- 
contamination detected between the vaccines. The tracheal 
swabs from Ark/B, Ark/C, Mass/B, Mass/C, and MassD/C vac- 
cinated birds were positive for IBV at all sampling times. 
The tracheal swabs from Del/C, GA/C, GA/A, and Mass/A 
vaccinated birds were positive on days 3, 6, and 9dpv but 
not at necropsy (14 dpv Expt. 1 or 13 dpv Expt. 2). Tracheal 
swabs for all of the contact-exposed birds in groups Ark/B, 
Conn/B, Mass/A, and Mass/B were positive for IBV at all sam- 
pling times. Swabs from Del/C, GA/A, Mass/C, and MassD/C 
contact-exposed birds were positive for IBV by RT-PCR or 
real time RT-PCR except at 3 dpv. Tracheal swabs from Ark/C 
contact-exposed birds were negative at 3, 6, and 9dpv but 
positive at necropsy (14 dpv Expt. 1 or 13 dpv Expt. 2). Using 
RT-PCR analysis, tracheal swabs collected from Ark/A vac- 
cinated and contact-exposed birds were negative for IBV. 
However, real-time RT-PCR analysis, which is more sensitive 
than RT-PCR, showed that the swabs from Ark/A vacci- 
nated birds but not contact-exposed birds were positive 
for virus. 


Table 3. Pairwise comparisons of amino acids between Conn vaccines from different manufacturers 
Residue position Conn/B vial 1 Conn/B vial 2 Conn/C vial 1 Conn/C vial 2 
10 Ala Thr Thr Thr 
58 Asn Asn Asn Deleted 
74 Ser Asn Asn Asn 
81 Thr Ser Thr Thr 
88 Asp Asp Gly Asp 
114 His His Tyr His 
161 Leu Phe Phe Phe 
188 Lys Asn Asn Asn 
214 Asp Asn Asp Asp 
243 lle Thr Thr Thr 
260 Val lle lle lle 
264 Phe Leu Leu Leu 
265 Thr Lys Lys Lys 
268 Tyr Asn Asn Asn 
270 Ser Thr Thr Thr 
281 Pro Leu Leu Leu 
309 Ser Gly Gly Gly 
330 Leu Pro Pro Pro 
331 Glu Lys Glu Glu 
343 Ser Leu Ser Ser 
3/3 Glu Gly Gly Gly 
376 Leu Ser Ser Ser 
381 Val Val lle Val 
383 Ser Leu Leu Leu 
387 Asp Lys Lys Lys 
388 His Ser Ser Ser 
389 Asn Asp Asp Asp 
402 Gly Asp Asp Asp 
406 lle lle Thr Thr 
430 Asp Asp Gly Asp 
460 Leu Met Met Met 
478 Tyr Tyr Cys Tyr 
481 Asn Thr Thr Thr 
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Table 4 Pairwise comparisons of amino acids between Mass vaccines from different manufacturers 


Residue position Mass/A Mass/B Mass/C MassD/C vial 1 MassD/C vial 2 
19 Ala Val Val Ala Ala 
38 Asp Asn Asn Asp Asn 
64 Gly Glu Gly Gly Gly 
69 Ile Thr Thr lle Thr 
73 Gly Asp Gly Gly Gly 
75 Val Val Val Val Ala 
76 Val Val Val Val Ala 

101 Tyr His His Tyr His 

109 Val Val Val Val Glu 

130 Ser Phe Phe Ser Phe 

179 Ala Glu Ala Glu Glu 

205 Arg Lys Lys Arg Lys 

246 Thr Ile lle Thr lle 

271 Asn Tyr Tyr Asn Tyr 

273 Thr Ser Ser Thr Arg 

376 Gly Glu Glu Gly Glu 

479 Ser Gly Gly Ser Gly 


Table 5 Pairwise comparisons of amino acids between GA 


vaccines from different manufacturers 


Residue position GA/A GA/C 
3 Gly Val 
52 Arg Gly 
58 Asp Ser 
60 Glu Asp 
68 Gly Asp 
84 Gly Arg 
86 Arg Ser 
90 Gln Lys 
a y2. Thr Ser 
22; Lys Gln 
141 lle Arg 
144 Asn Asp 
163 Asn Ser 
197 Val Ser 
198 Met Tyr 
199 Arg Ala 
200 Asp Arg 
201 Thr Ser 
257 His Asp 
282 Glu Asn 
283 His Gln 
330 Gly Arg 
389 Asp Glu 
390 Val Phe 
391 Asn Asp 
409 Arg lle 
423 lle Deletion 
479 Lys Asn 
481 Thr lle 


The Conn/B vaccine virus was only reisolated from 
the swab taken at 3dpv, from vaccinated birds while the 
virus was reisolated from all four sampling times from the 
contact-exposed birds (Table 6). From the birds that were 
vaccinated with Conn/C vaccine, which also contains the 
Mass virus, we detected IBV by RT-PCR, but Conn specific 
primers failed to amplify the Conn virus indicating that 
the viruses detected in the vaccinated and contact-exposed 
birds were not Conn virus (data not shown). Sequence anal- 
ysis confirmed that only the Mass virus was reisolated. 


S1 sequence analysis of reisolated vaccine viruses 


The $1 amino acid sequences from the reisolated viruses 
were compared to the sequence from the same vaccine virus 
that was used for vaccination. Analysis of the virus reiso- 
lated from Del/C, GA/C, Mass/A and Mass/C vaccinated and 
contact-exposed birds showed that the $1 gene amino acid 
sequence did not change after the virus replicated in birds. 

The S1 amino acid sequence analysis of vaccine viruses 
reisolated from Ark vaccinated and contact birds showed a 
consistent amino acid change of Tyr326 to Asn326 and a dele- 
tion of Asn345 (Table 7) in the vaccine virus obtained from 
all three manufacturers. Analysis of the original vaccine 
sequences revealed minor nucleotide peaks in otherwise 
clean sequence at Tyr326Asn, as shown in Fig. 1 as well as 
for Tyr43His and Ser213Ala, indicating multiple RNA species 
in the Ark/A and Ark/C vaccines. The deletion at amino acid 
residue 345 was not detected in the original vaccines. In 
addition, amino acid changes Ser130Gly (Ark/A and Ark/C) 
and Thr323Arg (Ark/C) detected in the reisolated vaccine 
viruses appeared to result from in vivo mutations since minor 
peaks were not detected in the original vaccine sequences, 
although selection cannot be excluded. 

The S1 amino acid sequence for the virus reisolated 
from Conn/B vaccinated birds had no changes when com- 
pared to the virus in the vaccine vial. The Conn/B virus 
obtained from contact-exposed birds had an amino acid 
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Table 6 RT-PCRand real time RT-PCR analysis of vaccine viruses reisolated from the trachea of vaccinated and contact-exposed 
birds at 3, 6, 9, 13 or 14 days post-vaccination 


Vaccines? 3 dpv° Tracheal swab 6 dpv Tracheal swab 9 dpv Tracheal swab 13 or 14dpv Tracheal swab 
Vaccinated Contact Vaccinated Contact Vaccinated Contact Vaccinated Contact 

Ark/A =e = —/+ = —/+ = —/+ = 
Ark/B + —/+ + + + + —/+ + 

Ark/C —/+ -_ + - —/+ — —/+ —/+ 
Conn/B + —/+ - + - + - + 
Conn/C — - - — - - - - 

Del/C + = + + + + = alls 
GA/A + = —/+ —/+ + —/+ — —/+ 
GA/C + + —/+ + —/+ + - - 
Mass/A —/+ + + + + — —/+ 
Mass/B + —/+ a + + —/+ + 
Mass/C + - + —/+ + + —/+ —/+ 
MassD/C + = + dp + tp + + 


4 All samples from negative groups were negative. 

5 dpv=days post-vaccination. 

© (—/+) Samples negative by RT-PCR/and positive by real time RT-PCR. 
d (—) Samples negative by RT-PCR and by real time RT-PCR. 


Table 7 Sequence of reisolated vaccine viruses that showed S1 gene sequence differences from the original vaccine sequence 


Vaccine Nucleotide position Nucleotide change Amino acid change Amino acid deletion 


Ark/A? 127% TSC Tyr43His 
388 A>G Ser130Gly 
637 T=G Ser213Ala 
976 ToA Tyr326Asn 
1033-1035 AAT deletion Asn345 
Ark/B® 9761+ TOA Tyr326Asn 
1033-1035! AAT deletion Asn345 
Ark/CS 127 TSC Tyr43His 
388 A>G Ser130Gly 
637! T>G Ser213Ala 
968! C—+G Thr323Arg 
976 THA Tyr326Asn 
1033-1035 AAT deletion Asn345 
Conn/B¢ 172—174+ AAT deletion Asn58 
238 TA No 
241 TSA No 
1033 TSC Ser345Pro 
GA/AS 662 GoT Ser221Thr 
988 GA Gly330Arg 
1152 CoT No 
1438 GT Lys479Asn 
Mass/Bf 1911+ A>G Glu6é4Gly 
313? TSC Ser105Pro 
1469 CT Pro490Leu 
4 Virus isolated at 13 days post-vaccination (dpv). 
> Virus isolated from vaccinated birds at 9 dpv and from contact birds at 13 dpv. 
¢ Virus isolated at 6 dpv. 
4 Virus isolated from vaccinated birds at 3 dpv and from contact birds at 6 and 13 dpv. 
© Virus isolated at 3 and 9 dpv. 
f 


Virus isolated from vaccinated birds at 3 and 9 dpv and from contact birds at 6 and 13 dpv. 
¥ Vaccinated birds. 
Contact birds. 
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Ark nucleotide 976 GA nucleotide 988 
Minor peak = A Minor peak= A 
A + 
GTTIATTAT TTTTGGGCT 


yf 
iil A 
iat 


} 


Tyr326Asn 


| 
1] 
| 


Gly330Arg 


Figure 1 


GA nucleotide 662 Mass nucleotide 64 
No minor peak Minor peak= G 


v v 


GACAGTCC CTGAGTG 


Ser221 Thr Glu64Gly 


Chromatograms of vaccine viruses showing major and minor peaks (Tyr326Asn, Gly330Arg, and Glu64Gly) indicating 


subpopulation selection associated with the reisolated vaccine viruses, and a single peak (Ser221Thr) indicating a mutation was 


associated with the reisolated vaccine virus. 


change of Ser345Pro, a deletion of Asn58 as well as two 
synonymous nucleotide changes (Table 7). Only one change, 
Ser345Pro in Conn/B appeared to result from in vivo muta- 
tions since no minor peaks were observed in the sequence 
chromatogram. 

The GA/A vaccinated birds had four nucleotide changes 
that resulted in three amino acid changes, Ser221Thr as 
shown in Fig. 1, as well as Gly330Arg and Lys479Asn. The 


Table 8 


Reisolated vaccine Amino acid change 


Gly330Arg change resulted from selection of a subpopula- 
tion since a minor A peak was observed in the chromatogram 
(Fig. 1), whereas the other changes appear to be the result 
of in vivo mutations. We also observed a synonymous C to T 
nucleotide change at position 1152. 

The S1 sequence in the virus obtained from Mass/B vac- 
cinated and contact birds had a Glu64Gly change. A minor G 
peak was observed in the chromatogram indicating that the 


Subpopulations and changes following passage of reisolated virus in embryonating eggs 


Double peaks? Egg passage change? 


Ark/A (vaccinated 13 dpv)‘° Tyr43His 

Ser130Gly 
Ser213Ala 
Tyr326Asn 


Asn345 deletion 


Tyr326Asn 
Asn345 deletion 


Tyr326Asn 
Asn345 deletion 


Tyr43His 
Ser130Gly 
Ser213Gly 
Thr323Arg 
Tyr326Asn 
Asn345 deletion 


Ark/B (contact 13 dpv) 


Ark/B (vaccinated 9 dpv) 


Ark/C (vaccinated 6 dpv) 


Asn58 deletion 
Ser345Pro 


Ser221Thr 
Gly330Arg 
Lys479Asn 


Glu64Gly 


Glu64Gly 
Ser105Pro 
Pro490Leu 


Conn/B (contact 14 dpv) 


GA/A (vaccinated 9 dpv) 


Mass/B (vaccinated 3 dpv) 
Mass/B (contact 13 dpv) 


Yes No 

No No 

Yes No 

Yes No 

- No 

No Yes (EP3)¢ 
— Yes (EP3) 
No No 

- No 

Yes No 

No No 

Yes No 

No No 

Yes No 

— No 

— No 

No Yes (EP4) 
No No (EP3) 
Yes No (EP3) 
No No (EP3) 
Yes No 

Yes No 

No Yes (EP3) 
No Yes (EP3) 


4 Major and minor chromatogram peaks observed in clean sequence data. 

> Viruses were passaged 10 times in 9—11-day old embryonating eggs unless otherwise noted. 

© Virus reisolated from vaccinated or contact-exposed birds on the day indicated (dpv = day post-vaccination). 
4 Each egg passage was sequenced. (EP#) = Egg passage number where the change was observed. 
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change resulted from selection of a subpopulation (Fig. 1). 
The virus isolated at necropsy from the contact-exposed 
birds showed additional changes of amino acid Ser105Pro 
and Pro490Leu (Table 7). Those changes apparently resulted 
from in vivo mutations. 

The S1 amino acid sequence analysis of MassD/C vacci- 
nated and contact-exposed birds were not consistent. The 
virus obtained from the first swabs from the vaccinated birds 
at 3 dpv and the contact-exposed birds at 6dpv had amino 
acid changes that were not observed in the viruses isolated 
from subsequent swabs (data not shown). 


Sequence analysis of reisolated vaccines passaged 
in eggs 


To determine if genetic changes observed in the reiso- 
lated vaccine viruses would revert to the original vaccine 
sequence we passaged them a maximum of ten times in 
embryonating eggs and examined selected passages for 
genetic reversion (Table 8). No reversion was observed in the 
Ark/A, Ark/B, and Ark/C viruses reisolated from vaccinated 
birds following 10 passages. The Ark/B viruses reisolated 
from contact-exposed birds reverted from Asn to Tyr at posi- 
tion 326 and an additional Asn at position 345 was observed 
at pass 3 in embryonating eggs, which was the same as the 
original vaccine. The Conn/B virus reisolated from contact- 
exposed birds reverted to the original vaccine sequence at 
position Pro345Ser after 4 passages in eggs but the dele- 
tion of Asp58 in the reisolated vaccine was still present 
following 10 passages. None of the GA/A reisolated vac- 
cine viruses reverted following 3 egg passages. The GA/A 
virus was only passaged 3 times in eggs due to contamina- 
tion. The Mass/B detected by RT-PCR from vaccinated birds 
at 9dpv could not be rescued in embryonating eggs; thus, 
the 3dpv swabs were passaged in eggs. Reisolated Mass/B 
from vaccinated and contact-exposed birds did not revert 
to the original vaccine sequence at position Glu64Gly fol- 
lowing 10 passages in eggs. Mass/B virus reisolated at 13 dpv 
from contact-exposed birds reverted to the original vaccine 
sequence at position Pro105Ser and Leu490Pro following 3 
passages in eggs. 


Discussion 


In this study, we examined the S1 sequence of 12 vac- 
cines for IBV and compared that sequence data to the same 
vaccines reisolated from vaccinated and contact-exposed 
chickens. We focused on the S1 gene for this analysis because 
it is widely known that the S1 subunit of spike contains 
serotype specific epitopes, as well as epitopes that induce 
virus-neutralizing antibodies and S1 is directly involved in 
attachment to host cells. In addition, it was shown that 
by far, most genetic changes occur in the $1 gene during 
adaptation to the host [15—17]. Until now, it was not clear 
if mutations or selection of a more fit subpopulation was 
responsible for the changes observed when coronaviruses 
are attenuated or adapted to a particular host system [16]. 
Our data clearly shows that selection of a subpopulation of 
intra-vaccine quasispecies results from infection and repli- 
cation of IBV vaccine viruses in chickens and that genetic 
mutation among IBV vaccines also occurs. 


With a couple of exceptions where it was necessary to 
sequence cloned amplicons due to sequencing difficulties, 
we sequenced RT-PCR products because that data repre- 
sents all of the genetic diversity within the population of 
viruses in a given isolate. Sequence ambiguities in otherwise 
clean sequence data can be used to identify diverse viral 
RNA genomes within an isolate [5]. Furthermore, the area 
under the peaks when two or more peaks are present at any 
single nucleotide position in the sequence chromatogram 
can be used to identify the major and minor subpopulations. 
Our data show differences in $1 sequence within the same 
vaccine serotype from different manufacturers, differences 
in S1 sequence between different vaccine serials from 
the same manufacturer, and intra-vaccine $1 differences. 
Differences in the S1 gene between the same vaccine 
serotype from different manufacturers likely reflects 
the source of the original virus strain used to produce 
the vaccine and the methods and number of passages in 
embryonating chicken eggs used to attenuate the virus. 
S1 gene differences between vaccine serials from the 
same manufacturer are most likely due to the generation 
of mutations through replication of the viral genome. 
The World Organization for Animal Health (OIE) guidelines 
(http: //www.oie.int/eng/normes/mmanual/A_summry.htm) 
for production of IBV vaccines stipulates that no more than 
5 passages beyond the master seed can be used to pro- 
duce the vaccine. IBV has been reported to have a 1.5% 
mutation rate per year, and SARS-CoV was reported to have 
a non-synonymous substitution rate of 1.16—3.30 x 10-3 
per site per year [7,8]. In addition, it has been reported 
that IBV can adapt to cell culture following as few as 2—5 
passages and as many as 25 amino acid substitutions have 
been documented in the S1 gene after only 7 passages 
in cell culture [18,19]. Clearly, coronaviruses have the 
capacity for rapid genetic change, which likely accounts for 
sequence differences between different serials of the same 
vaccine. 

To our knowledge, this is the first report where both 
selection of a more fit virus subpopulation and mutations 
are associated with replication of modified live IBV vaccine 
viruses in chickens. Vaccine virus subpopulation selection 
following vaccination, as evidenced by major and minor 
peaks in otherwise clean sequence data, were observed 
for Tyr43His, Ser213Ala, and Tyr326Asn changes in reiso- 
lated Ark/A and Ark/C vaccine, Gly330Arg in reisolated 
GA/A vaccine and Glu64Gly in reisolated Mass/B vac- 
cine. Genetic mutations as evidenced by single peaks in 
clean sequence data and no reversion following egg pas- 
sage of reisolated viruses, were observed for Ark/A at 
Ser130Gly, and Ark/C at Ser130Gly and Thr323Arg. Muta- 
tions were also detected for GA/A at positions Ser221Thr 
and Lys479Asn. However, it should be noted that exist- 
ing viral RNA species could be responsible for apparent 
genetic mutations if that viral RNA subpopulation was below 
our level of detection. The level of detection for the 
RT-PCR test used herein is approximately 500 template 
copies/ml (unpublished data). No genetic subpopulations 
for Conn/B at position Ser345Pro were detected in the 
original vaccine, however; the reisolated vaccine reverted 
to Ser after only 4 egg passages. This was also observed 
for Mass/B at positions Ser105Pro and Pro490Leu, which 
reverted after only 3 egg passages. Reversion after only 
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a few egg passages suggests that the genetic change is 
important for growth in eggs but it is not clear if selection 
of a subpopulation or mutations was responsible for those 
changes. 

A deletion of 3 nucleotides (AAT) coding for Asn345 was 
observed in all of the reisolated Ark vaccines indicating 
that that genetic change is important for in vivo repli- 
cation of Ark type vaccines. However, based on sequence 
data from cloned RT-PCR product, the reisolated Ark/B vac- 
cine reverted following embryonating egg passage indicating 
that a low-level virus population retaining the 3 nucleotide 
sequence was being maintained in vivo, since it is unlikely, 
but not impossible, that the virus could acquire an AAT inser- 
tion through a recombination event or other mechanism 
after only a few back passages in eggs. It was necessary to 
clone the RT-PCR products from the Ark/B vaccine to obtain 
sequence data of single fragments, presumably because 
viruses with and without the deletion were being maintained 
in the population. 

Amino acid sequence differences in the IBV spike glyco- 
protein are responsible for diverse serotypes of the virus, 
which do not cross-protect [20]. Some IBV vaccines, so called 
protectotypes, have been shown to provide more cross- 
protection than other similar vaccines [21]. It is likely that 
these protector-type vaccines have a more diverse popu- 
lation of viruses capable of simulating a broader immune 
response in the host. In addition, it has been shown that viral 
subpopulations can work together to invade and replicate 
in the host [22,23]. Thus, it seems logical that a vac- 
cine having a more dynamic quasispecies population ought 
to be more efficacious than one producing fewer mutant 
sequences. However, molecular evolution can lead to the 
selection of virulent viruses as well as the emergence of 
new viral pathogens [9]. A largely undocumented but widely 
accepted phenomenon of rolling vaccine reactions in com- 
mercial chicken flocks occurs when less than 100% of the 
flock receives the vaccine allowing repeated transmission 
of the vaccine virus from vaccinated birds to unvaccinated 
flock-mates. The highly infectious nature of IBV combined 
with vaccination and re-infection of vaccine viruses in a 
commercial chicken flock that can contain over 20,000 birds, 
likely allows the virus to undergo molecular evolution. The 
evidence for this can be found in the multitude of variant 
viruses with sequences extremely similar to the vaccines 
used in the flocks [24]. 

In summary, it appears that high-titered IBV vaccines 
contain a mixture of genetic variants and vaccination of 
chickens results in selection of the most fit viral RNA subpop- 
ulations as well as generation of mutations. This allows rapid 
adaptation to the host ensuring replication and stimulation 
of an immune response, however; accumulating evidence 
also indicates that it can lead to persistence and the emer- 
gence of new viral pathogens [5,9]. Recognizing selection 
and mutation of vaccine viruses in the host is important 
for our understanding of how vaccine viruses can poten- 
tially replicate, persist and evolve to cause disease, and 
applies to all modified live coronavirus vaccines including 
potential vaccines against the SARS-CoV. This information 
can also be used to develop safer, higher quality, more 
targeted modified live vaccines containing viral subpop- 
ulations that successfully induce an efficacious immune 
response. 
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